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Abstract. Using a ladder-rainbow kernel previously established for light quark hadron physics, 
we explore the extension to masses and electroweak decay constants of ground state pseudoscalar 
and vector quarkonia and heavy-light mesons in the c- and b-quark regions. We make a systematic 
study of the effectiveness of a constituent mass concept as a replacement for a heavy quark dressed 
propagator for such states. The difference between vector and axial vector current correlators is 
explored within the same model to provide an estimate of the four quark chiral condensate and the 
leading distance scale for the onset of non-perturbative phenomena in QCD. 
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DYSON-SCHWINGER EQUATIONS OF QCD 

A great deal of progress in the QCD modeling of hadron physics has been achieved 
through the use of the ladder-rainbow truncation of the Dyson-Schwinger equations 
(DSEs). The DSEs are the equations of motion of a quantum field theory. They form 
an infinite hierarchy of coupled integral equations for the Green's functions (n-point 
functions) of the theory. Bound states (mesons, baryons) appear as poles in the appro- 
priate Green's functions, and, e.g., the Bethe-Salpeter bound state equation appears after 
taking residues in the DSE for the appropriate color singlet vertex. For recent reviews 
on the DSEs and their use in hadron physics, see Refs. [m,a,BL^. 

In the Euclidean metric that we use throughout, the DSE for the dressed quark 
propagator is 

Sip)-^ = Z2i/p + Z^m{ll)+Zi j\^D^y{p-q)^Y^S{q)r^{q,p) , (1) 

where D^v{k) is the renormalized dressed-gluon propagator, ry{q,p) is the renormal- 
ized dressed quark-gluon vertex. The solution of Eq. ([T]) is renormalized according to 
"^(p) ^ = ^7' P + f^il^) ^ sufficiently large spacelike /i^, with m(/i) the renormalized 
quark mass at the scale /i. We use jl = 19GeV for numerical work. The renormalization 
constants Z2 and Z4 depend on the renormalization point and the regularization mass- 
scale. 

After taking the residue at the bound state pole of the inhomogeneous DSE for the 
relevant vertex, bound state Bethe-Salpeter equation is 

r'Hp+.P-) = rK{p,q;P)S%q+)rHq+,q-)S\q-) , (2) 
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FIGURE 1. The effective ladder-rainbow kernel from the Maris-Tandy [MT] model [5] showing the 
infrared (IR) and ultraviolet (UV) components. 



where K is the renormalized qq scattering kernel that is irreducible with respect to a pair 
of qq lines. The quark momenta are ^+ = ^ + r\P and q- = q— {\ — r\)P where r\ is the 
momentum partitioning parameter. The choice of rj is equivalent to a choice of relative 
momentum q\ physical observables should not depend on r\ . This provides us with a 
convenient check on numerical methods. The meson momentum satisfies = —m^. 



RAINBOW-LADDER TRUNCATION 

A viable truncation of the infinite set of DSEs should respect relevant (global) symme- 
tries of QCD such as chiral symmetry, Lorentz invariance, and renormalization group 
invariance. For electromagnetic interactions and Goldstone bosons we also need to re- 
spect color singlet vector and axial vector current conservation. The rainbow-ladder 
truncation, which achieves these ends, replaces the BSE kernel by the (effective) one- 

gluon exchange K{p,q;P) — > -4n aef{{k^) D^^^{k)YY^ ® y7v along with the replace- 
ment of the DSE kernel for S{p) by Zig^D^y{k)ry{q,p) ^ 47caeff{k^)D^/[^^{k)Yv^ 
where k = p — q, and OCeff (^^) is an effective running coupling. This truncation is the 
first term in a systematic expansion [0, M\ of the quark-antiquark scattering kernel K; 
asymptotically, it reduces to leading-order perturbation theory. Furthermore, these two 
truncations are mutually consistent in the sense that the combination produces color sin- 
glet vector and axial-vector vertices satisfying their respective Ward identities. In the 
axial case, this ensures that the chiral limit ground state pseudoscalar bound states are 
the massless Goldstone bosons associated with chiral symmetry breaking [9, 10]. In the 
vector case, this ensures, in combination with impulse approximation, electromagnetic 



current conservation MM- 

We employ the ladder-rainbow kernel found to be successful in earlier work for light 
quarks LIQ,^. It can be written aeff(^^) = oc^^ik^) + a^^{k^) The IR term implements 



TABLE 1. DSE results |5] for pseudoscalar and vector meson masses and electroweak decay 
constants, together with experimental data |0]. Units are GeV except where indicated. Quantities 
marked by t are fitted with the indicated current quark masses and the infrared strength parameter 
of the ladder-rainbow kernel. 
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calc 


0.138 0.131 
0.138''' 0.131''' 


0.496 0.160 
O.497I' 0.155 


0.770 0.216 
0.742 0.207 


0.892 0.225 
0.936 0.241 


1.020 0.236 
1.072 0.259 



the strong infrared enhancement in the region < < 1 GeV^ required for sufficient 
dynamical chiral symmetry breaking. The UV term preserves the one-loop renormaliza- 
tion group behavior of QCD: aeff(^^) as{k^y^°°P (k^) in the ultraviolet with Nf = 4 
and Aqcd = 0.234 GeV. The strength of a}^ along with two quark masses are fitted to 
(qq), mjj^j]^ and /;;:. This and selected light quark vector meson results are displayed in 
Table[T] The momentum dependence of the kernel along with the IR and UV components 
is displayed in Fig.[TJ The infrared component of this effective kernel is phenomenolog- 
ical because QCD is unsolved in such a non-perturbative domain. To help replace such 
phenomenology by specific mechanisms, it is necessary to first characterize its perfor- 
mance in new domains. 

RESULTS FOR HEAVY QUARK MESONS 

In Table [2] we display the results for the heavy-light ground state pseudoscalars and 
vectors involving a c-quark or b-quark. We use DSE solutions for the dressed light 
quarks. If a constituent mass propagator is used for the heavy quark, with the constituent 
mass obtained from a fit to the lightest pseudoscalar, the various meson masses are 
easily reproduced. The constituent masses found this way are M^°"^ = 2.0 GeV for 
the c-quark, and M^°"^ = 5.3 GeV for the b-quark. To compare with what is known 
about quark masses, we take the quark current mass values [6] = 1.2 ±0.2 GeV, and 
nit = 4.2 ± 0.2 GeV at scale /i = 2 GeV and use the quark DSE to run the masses into 
the timelike region where the meson mass shells are located. If all the meson momentum 
runs through the heavy quark, the mass function from the DSE would suggest an 
effective quark mass Mq^^{p^ ~ —M^) where M is the meson mass. In this domain 
the mass function varies slowly and reproduces the previously obtained values for M^°^^ 
within 10%. In this sense, heavy quark dressing is well summarized by a constituent 
mass. However, the electroweak decay constants obtained from the constituent mass 
approximation are 30-50% below the available experimental values. Moreover, within 
this ladder-rainbow model, quark dressing is not a minor effect because the use of fully 
dressed quark propagators, both heavy and light, does not yield a physical bound state 
solution for these heavy-light states involving a c-quark or b-quark. 

The results for equal quark mesons (quarkonia), displayed in Tabled show a different 



TABLE 2. Calculated masses and electroweak decay constants for ground state pseudoscalar 
and vector heavy-light mesons, together with experimental data [6], all in GeV. In the rows 
labelled calc M, the heavy quark is described by a constituent mass fit to the lightest pseudoscalar 
(marked by t)- In the rows labelled Euv only, the heavy quark is dressed, but only by the UV 
component of the DSE kernel; this did not produce physical solutions for D or D* states. 





D D* 


D, D* 


B B* 


B, B* 


B, B* 


expt M 
calc M 
Euv only 


1.86 2.01 
1.85^ 2.04 


1.97 2.11 
1.97 2.17 


5.28 5.33 
5.27'!' 5.32 
4.66 


5.37 5.41 

5.38 5.42 
4.75 


6.29 ? 
6.36 6.44 
5.83 


expt f 
calc f 
Euv only 


0.222 ? 
0.154 0.160 


0.294 ? 
0.197 0.180 


0.176 ? 
0.105 0.182 
0.133 


7 7 
0.144 0.20 
0.164 


7 7 
0.210 0.18 
0.453 



perspective. With the same fitted constituent masses, the meson masses are well repro- 
duced, and again the electroweak decay constants are too low by some 30-70 %. Use 
of dynamically dressed propagators removes almost all of this deficiency and the decay 
constants are well reproduced. This improvement provided by dynamical dressing of c- 
and b-quarks is persistent and systematic in the following sense. When the dressing is 
progressively introduced into all three stages (bound state solution, normalization loop 
integral, and then the loop integral for evaluation of the decay constant), the result al- 
ways increases towards the experimental value. This suggests that a constituent mass 
approximation, even for b-quarks, is inadequate. Small departures from a strictly con- 
stant mass function and field renormalization function Z{p^) for quarks in the relevant 
region of the complex plane are magnified due to the very weak binding of the mesons 
in question. 

With increasing quark mass, the quarkonia states become smaller in size and the 
ultraviolet sector of the ladder-rainbow kernel becomes more influential. The size of 
the dressed quark quasi-particle also decreases with quark mass so that the entire heavy 
quarkonia dynamics should be dominated by the ultraviolet sector of the kernel as 
dictated by pQCD. In the last two rows of Tabled we display the results obtained by 
using just the infrared sector, or just the ultraviolet sector, of the ladder-rainbow kernel 
in all phases of the calculation, binding and quark dressing. For ease of calculation, we 
define these two components of the kernel as shown in Fig. [U corresponding to the two 
terms in the kernel expression [[loL [sj]. With only the UV kernel, a mass shell couldn't 
be reached for cc, but excellent masses for the bb states are obtained. 

In the case of heavy-light mesons, the binding and dressing effect for the light quark is 
not significantly different from the dynamics of light quarkonia and the empirical chiral 
condensate that set the characteristic infrared length scale and strength of the employed 
ladder-rainbow kernel. For this reason we investigated a hybrid procedure in which the 
ultraviolet component of the kernel was used for heavy quark dressing, while the full 
kernel was used for binding and light quark dressing. The results displayed in Tabled 
show that physical B states are produced this way, whereas that was impossible with a 
b-quark dressed via the full kernel. This again suggests that for heavier mass, or smaller 
size objects, the infrared sector of the present kernel is too strong. In the case of heavy 
quarkonia, the bound state wavefunction itself provides some limitation of support in the 



TABLE 3. Calculated masses and electroweak decay constants for ground state pseu- 
doscalar and vector quarkonia, together with experimental data |6], all in GeV. Results em- 
ploy either a constituent mass treatment of the heavy quark, labelled Mg°"^ or a dynamical 

dressed propagator, labelled Eg^^(/:i^). The last two rows shown the % change from the full 
dynamical calculation due to retention of just the IR or the UV component of the ladder- 
rainbow kernel in all aspects of the calculation, both quark dressing and binding. 







^y/v/ fj/v 


Mm, fnt 


Mt /t 


expt 

calc with Mg'"' 
calc withEg^^^(/?2) 


2.98 0.340 
3.02 0.239 
3.04 0.387 


3.09 0.411 
3.19 0.198 
3.24 0.415 


9.4 ? ? 
9.6 0.244 
9.59 0.692 


9.46 0.708 

9.65 0.210 

9.66 0.682 


IR only 
UV only 


-21% -29% 


-20% -20% 


-15% -75% 
-0.6% -12% 


-16% -50% 
-0.8% -25% 



BSE integrand to small distance. However that is less the case for the heavy-light states; 
a forced suppression of the infrared sector of the heavy quark DSE seems necessary. 

This may indicate the present kernel, determined for light quarks, has implicit infrared 
strength through quark-gluon vertex dressing that should be much reduced for heavy 
quarks. Aspects of this study may also be a partial confirmation of Brodsky and Shrock's 
suggestion of a universal maximum wavelength of quarks and gluons in hadrons [|l2ll . 
The resulting minimum momentum would increase with mass of the hadron or dressed 
quark/gluon constituent and such a cutoff may systematically link our findings. 



LEADING NPQCD SCALE AND THE FOUR QUARK 

CONDENSATE 

Quark helicity and chirality in QCD are increasingly good quantum numbers at short 
distances or at momentum scales significantly larger than any mass scale. One manifes- 
tation of this is that, for chiral quarks, the correlator of a pair of vector currents is iden- 
tical to the corresponding correlator of a pair of axial vector currents to all finite orders 
of pQCD. In the non-perturbative circumstance, the difference of such correlators mea- 
sures chirality flips, and the leading non-zero contribution in the ultraviolet identifies the 
leading non-perturbative phenomenon in QCD. This is the four quark condensate [fH]. 
The ladder-rainbow kernel can produce the difference correlator as a vacuum polariza- 
tion integral in momentum space, where the propagators are dressed and the vector and 
axial vector vertices are generated in a way consistent with the symmetries. We use the 
large spacelike momentum dependence to extract the leading coefficient or condensate; 
and a Fourier transform identifies the leading non-perturbative distance scale. 
The vector current-current correlator is formulated as the loop integral 

Til,{P) = j\'xe'-{0\Tj^{x)j+{Om = - j -^Jr{y^S{q+)Yl{q,P)S{q^)] , 

(3) 

where A indicates regularization, e.g., by the Pauli-Villars method, and Ty is 
the dressed vector vertex. The axial vector correlator is formulated in an anal- 




ogous way and we directly calculate the difference correlator which does not 
require ultraviolet regularization. With nj^y(/') = (P^S^y — PjiPv 

)n^(/'2), and 

I{^^{P) = {P'^5^v-PpiPv)Yl^{P'^)+P^PvYlj^{P'^), the quantity of interest here is 

nj-^(p2) =n^(p2)_nA(p2)_ 

The leading non-perturbative contribution to starts with dimension d = 6 and 



involves the four-quark condensate in the form [|14l Il5|] 



nr(^V-^^{.+2^lf + .n(^)l} + 0(i,) . (4, 

The left panel of Fig. [2] indicates that our numerical calculation of P^Y]^^^{P^) identi- 
fies a leading ultraviolet constant reasonably well. The four quark condensate {qqqq) ex- 
tracted via Eq. (Hj) is 65% greater than the common vacuum saturation assumption {qq)^ 
at the renormalization scale /i = 19 GeV used in this work. The low P^ limit provides a 
reasonable account of the first Weinberg sum rule iI^IIt!] : D^^^ (P^) \p2^Q = —f^, in 



the /;;: = 0.0924 GeV convention. We obtain fj^ = 0.0728 GeV that way. Our resuks are 
consistent with the second Weinberg sum rule [16] P'^n\.-^{P^)\p2^^ = 0. The Das- 
Guralnik-Mathur-Low- Young sum rule [[Isl] relates Jq dP^ P^Uj^'^iP^) to the electro- 



magnetic component of mj^± — m^o. We obtain 4.86 MeV for this mass difference in 
comparison with 4.43 ±0.03 from experiment. 

The right panel of Fig. [21 displays the x-dependence of the V-A correlator amplitude 
expressed as a ratio to the free V+A correlator amplitude. The latter diverges in the 
ultraviolet as \n{P^), which corresponds to x^^ at small x. This chirality-flip ratio 
identifies a scale of ~ 0.5 fm for the onset of non-perturbative dynamics. This ladder- 
rainbow truncation within a DSE format produces the same scale as obtained from the 
ratio of P-S and F+S correlators in both a lattice-QCD calculation and the Instanton 
Liquid Model IH. 
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